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We analyze the initial value problem for spinor fields obeying the Dirac equa- 
tion, with particular attention to the characteristic surfaces. The standard Cauchy 
initial value problem for first order differential equations is to construct a solution 
function in a neighborhood of space and time from the values of the function on a 
selected initial value surface. On the characteristic surfaces the solution function 
may be discontinuous, so the standard Cauchy construction breaks down. For the 
Dirac equation the characteristic surfaces arc null surfaces. An alternative version 
of the initial value problem may be formulated using null surfaces; the initial value 
data needed differs from that of the standard Cauchy problem, and in the case we 
here discuss the values of separate components of the spinor function on an inter- 
secting pair of null surfaces comprise the necessary initial value data. We present 
an expression for the construction of a solution from null surface data; two ! ana- 
logues of the quantum mechanica 1 Hamiltonian operator determine the evolution 
of the system. 
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1. Introduction 

In reference [1] we discussed the Cauchy problem, initial values, and charac- 
teristics of the Klein Gordon equation for scalar fields. In this work wc extend 
the discussion to include spinor fields obeying the free particle Dirac equation. Al- 
though the approach is the same and frequent references are made to reference 1 
we have attempted to make this paper self-contained. 

The standard Cauchy problem for a first order partial differential equation is 
to use the values of the solution function on an initial value surface to determine 
the values of the function in a neighborhood of the surface [2]. Solution of the 
problem proceeds by showing how the initial value data and the differential equa- 
tion determine all the time derivatives of the function on the initial surface and 
thus allow a power series development of the function for future times. 

If the solution function has a discontinuity across some special surface then the 
standard Cauchy problem cannot be solved using that surface for the initial value 
data. Such surfaces are called characteristic surfaces or simply characteristics. The 
characteristics of the Dirac equation for free particles in fiat space are easily shown 
to be null surfaces. Thus the discontinuities of the solution functions propagate at 
the velocity of light, independently of the mass parameter in the equations. The 
characteristics of Maxwell's equations, the Klein Gordon equation, and the vacuum 
Einstein equations of general relativity are also all null surfaces [1], [3]. 

We first briefiy review the standard Cauchy problem for the Dirac equation. 
We then obtain the characteristics and solve the analogous problem for the pla- 
nar null surfaces. Initial data for this situation consists of the values of disparate 
components of the spinor function on the pair of null surfaces. We give a sim- 
ple evolution operator expression in terms of null coordinates, which is a direct 
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analogue of the usual expression of quantum theory, except that two analogues 
of the Hamiltonian operator determine the evolution of the solution, that is are 
generators of displacements in the null coordinates [4, 5, 6, 7, 8]. The result is 
remarkably similar to that obtained for the scalar field [1]. 

2. Development Of a Spinor Function 
From a Constant Time Hypersurface 

Since the Dirac equation may be expressed in Hamiltonian form the develop- 
ment of a solution from initial value data given on a constant time hypersurface is 
one of the most familiar problems of physics [9] . We express the Dirac equation in 
Hamiltonian form as [10] 

iip\t{x, t) = [-id. V +(3m]ip{x, t) = HiIj{x, t) (2.1) 

(We use units in which h = c = 1. The slash notation indicates differentiation 
with respect to the indicated variable, in this case time t.) We write a Taylor series 
expansion as 

ipix, t) = ip{x, t)t=o + tij^ti^, t)t=o + 2! ^)^=o + - (2-2) 

All of the time derivatives at t = may be readily obtained from the Dirac 
equation and the value of the spinor function at t = 0, which we call h{x). When 
substituted in the series these give the familiar result 
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( —iHt)'^ 

iPix, t) = i;{x, t)t=o+{-iHt)ilj{x, t)t=o+ ^ ^, ' ^(f, t)t=o+... = e-'"^h{x) (2.3) 

The indicated exponential of the Hamiltonian operator is thereby seen as the 
finite time displacement operator. 



3. Characteristics of the Dirac Equation 

The characteristics of the Dirac equation are found in the same way as for the 
Klein Gordon equation in [1], except that since the Dirac equation is first order in 
time the problem is even easier. We suppose that the initial value of the solution 
is given on a hypersurface S, t — T{x), by a four component spinor (see fig. 1) 

i;{x,T{x)) = h{x) (3.1) 
and that the solution is given by a Taylor series expansion 



i^ix, t) = i^ix, t)t=T +[t- T(f)]V'|t (f , t)t=T + i^^j^V^l^li (f , t)t=T + ... (3.2) 

To find the first time derivative in this we first express the gradient of the 
initial value h in terms of t = T{x), as 



V/i(f) = Vi^{x, t)t=T + iJ\t{x, t)t=T V T{x) (3.3) 

From this and the Dirac equation (2.1)we may write an equation for the time 
derivative of ^0 on S, 



i[I - a. V T{x)]'il;\t{x, t)t=T ioi. v Kx) + (3mh{x) (3.4) 

This may be solved for the time derivative if the matrix in square brackets has 
an inverse. The determinant of that matrix is easily found to be 

\I -a.VT{x)\^{l-VT{xff (3.5) 

Thus the condition that S be a characteristic is that the above quantity be 
zero, or 

Vr(f)2 = 1 (3.6) 

This is the same as the characteristic equation for the Klein Gordon equation; 
thus the characteristics are null surfaces, and are independent of the mass param- 
eter in the Dirac equation. The particular null surfaces that we will study in the 
remainder of this paper are 

u — t — X — Q u characteristic or null surface 

V — t -\- X — Q V characteristic or null surface (3.7) 
These are shown in figure 2. 
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4. Construction of a Solution Function 
in Terms of Null Coordinates 



The covariant form of the Dirac equation is [10] 

d 

T^^^V' = rni; (4.1) 

We will study this using the null coordinates (see fig. 2) u — t — x and 
V — t + X, and suppress dependence on y and z in all that follows [11]. Then the 
gamma matrices are ^'^ = 7^ — 7^ and 7^ = 7^ + 7^. The Dirac equation and the 
gamma matrix algebra in terms of the null coordinates are 



^ + ^ = 1, (7")' = (7^)' = (4.3) 



We define a pair of projection operators as 



K = ^ K^-^ (4.4) 



From (4.3)the usual projection operator properties follow, that is 

+ = 1, A„A„ = A„, A^A^ - A„, A^A^ = A„A„ = (4.5) 
To analyze the Dirac equation on and near the null surfaces u — Q and v — Q 
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we make an expansion in the mass m, since the solutions of the zero mass Dirac 
equation are functions of only u or only v. Thus we write 



,il;^^\u,v) + rrv 



(4.6) 



Substitution of this into the Dirac equation gives the following set of iterative 
equations 



Here the functions f and g are any differentiable 4-tuple functions of u and v; we 
will call these the generating functions. It is instructive to consider a representation 
of the gamma matrices in which the projection operators L are diagonal; then the 

contains I's at entry (1,1) and (2,2), and zero everywhere else; similarly, A^; has 
I's at entries (3,3) and (4,4) and zeros everywhere else. 

As for the zero order solution , they are (in transpose notation): 



Thus only two components of f and two components of g enter the solution. 



Solution of the zero order equation is easy by inspection. 



(4.8) 



(4.9) 
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The first order equation (4. 7) is 



[7"^|^ + 7"^|;]^^'H«,^) = Kf{u)+Av9{v) (4.10) 

Only a particular solution to this is needed since the homogeneous solution 
may be absorbed into the zero order solution. We break the equation into u and v 
parts by premultiplying by and using (4. 2) to get 

^^A^^|JW(u,v) = ^J'^A,g{v) (4.11) 

The right side of this is a function of v only, so the left side must be also. Thus 
we may write a solution by integration. 



Au^^'\u, v) = A. j dv'g{v') (4.12) 

Here vq is an arbitrary parameter. In the same way we may obtain a similar 
expression for the v projection of the solution, 

u 

K^^^\u, v) = ^7'A« j du'fiu') (4.13) 

Mo 

The sum of (21) and (22) gives the complete first order solution 

V u 

^W{u,v) = ^[7"A, / dv'g{v')+j^A^ [ du'f{u')] (4.14) 
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The second order equation (16) may now be written as 



.d_ 



-Aui^^^Hu, v) = ^[7"A, j dv'g{v') + 7"A„ j du'f{u')\ (4.15) 
As before we premultiply by to obtain, with use of (4.4)and (4.5), 

u 

|-A,V^(2) ^) ^ zIa^ y du'fiu') (4.16) 



Wo 



Since the right side of this is a function of only u we may write a particular 
solution by inspection, 

u 

Aui^^^\u, v)^^{v- vo)K J du'fiu') (4.17) 

Mo 

In the same way we may obtain the v projection of the solution, 

V 

Avi;^'^\u, v)^^{u- uo)Ay J dv'g{v') (4.18) 

vo 

Summing (4.17)and (4.18)we have the complete second order solution. 



^P^^\u,v) = ^{v-vo)[A^J du'f{u') + {u-uo)A,J dv'g{v')] (4.19) 

Wo Vo 



The procedure is now clear, and wc may continue to all orders; the even orders 
are similar in form, and the odd orders are similar in form. The complete series 
solution is 
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V'(m,u) = 2^ r«)A«/H+2^ ^ Tj (l-X'^ Ty)Kg[v) 

n=0 ■ j=0 

(4.20) 

where the multiple integral operator is defined as 

u u' m"-i 

r::/w = 1 j ... j du-fiu^) (4.21) 

Wo Wo Wo 

The series (29) may be readily summed to give a concise expression for the 
solution 

i;{u,v) = e-^(^-^°)r"(l - !^^-r„)A„/(«) + e-^(---°)r^(l - ^7"r,)A,^(^) 

(4.22) 

This gives a solution of the free particle Dirac equation for any pair of generat- 
ing functions f and g; the generating functions together have only four independent 
components however, so this initial value data contains the same amount of infor- 
mation as that in the standard Cauchy problem. 
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5. Initial Values on Null Surfaces 

The generating functions f and g are arbitrary functions, and are simply related 
to the initial values of the solution on the null surfaces as we will now discuss. In 
this and the following section we set uq — vq — without loss of generality. 

To get a relation between the initial values of the function y on the null surfaces 
and the generating functions f and g we set u = and then v — in (4.22)and 
easily find 

7777 

M^, 0) = (1 - — 7-r„)A„/(u) + r,A,5(o) (5.1) 

7171 

MO, ^) = (1 - ^7"r,) A,5(^) + r„A„/(0) (5.2) 

We wish to make a convenient choice for the values of /(O) and g{0); taking 
u — V — in (5.1), (5.2)we see that 

^0(0,0) =A«/(0)+A,(7(0) (5.3) 

Accordingly we choose /(O) = g{0) = ipo{0, 0); then the quantities that appear 
in the evolution expression (4.22)are 

tin 

(1 - — 7"r«)A«/(ii) = V'o(0, 0) - A^Vo(0, 0)) (5.4) 

Z7TI 

(1 - —j''T^)A^g{v) = 0) - AuMO, 0)) (5.5) 
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Thus the spinorial expressions in parentheses turn the arbitrary 4-tuple func- 
tions into the indicated initial data on the null surfaces. These expressions may be 
inverted if it is desired to obtain f and g. 

Kf{u) = (1 + — 7'T«)[^o(0,0) - A,^o(0,0)] (5.6) 

7777 

Kg{v) = (1 + _7«r,)[V'o(0,0) - A,Vo(0,0)] (5.7) 
In terms of the initial values of the solution the expression (31) now reads 

(5.8) 

This is a complete expression for the evolution of the solution from its values 
on the pair of null surfaces tt = and v = Note that only the square of the mass 
appears in the evolution operators. 

The form of the solution in (5.8)is very similar to that obtained in [1] for 
scalar fields obeying the Klein Gordon equation; the two exponential evolution 
operators are identical and only the forms of the initial data expressions are slightly 
different since the present one contains spin information. Thus in terms of the null 
coordinates the dynamics of scalar and spinor fields is remarkably similar, much 
more so than in terms of the usual Cartesian coordinates. 
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6. Plane Waves 



We wish to verify the consistency of (5.8)for plane wave solutions of the Dirac 
equation. We will show that if the appropriate initial data functions are put into 
(5.8)then the solution function generated is the appropriate plane wave. We write 
a plane wave solution in Cartesian and null coordinates as 

iPiu, v) = e-'^*^'^''w{E, k) = e-'^'^-'^''w{\ r) (6.1) 
Here w is the usual Dirac 4-tuple spin function and the null momenta are given 

by 

^ E + k E-k E^-e 2 ,«o^ 

^ — — 7^ — ) T = — 7^ — ) ^'^ = : = — = A* (6-2) 

2 2 4 4 ^ ^ 

Thus the appropriate initial value quantities that enter the evolution expression 
(5.8)are 

V^o(«, 0) - AvM^, 0) = (e-*^^ - Av)w{\, r) 

^^0(0, v) - AuMO, 0) = (e-^^^ - A^)w{X, r) (6.3) 

We substitute these expressions into (5.8), abbreviating for convenience B — 
C = -IT, -^^BC, and find 
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= (e^^^r^gBu ^ e^^^^''e^'')w(A, r) - Fc{uv)w{X, r) (6.4) 

Here denotes the function related to a constant generator, defined and 
evaluated as 

Fc(uv) = ±^-=i^ = e-4^'-l = e-^''''l (6.5) 



n=0 



Notice that this function is symmetric in u and v. The quantity in parentheses 
in (6.4)is straight forward to evaluate; we expand the exponentials in double series 

as 



^ ^ ^ > 2^ n\ j\ ^ m\ k\ 

n=Oj=0 m=0,k=0 

(6.6) 

It is easy to obtain the operation of F"^ on powers of u, 



^ (i + l)(j + 2)...(j + n) ^^-^^ 



We substitute this into (6.6)and rearrange summation indices to find 
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n<j m>k 



n=0,j=0 ■' m=0,fc=0 n=0 ^ ' 

(6.8) 

We now combine (40) and (44) to get the complete solution 



V;(m, v) = e-^'^+^^w;(A, r) = e-*^"-*^'^w;(A, r) (6.9) 

which is the correct plane wave solution. We have thus checked the consistency 
of a known solution with our formalism. Moreover it follows that the evolution 
operator equation (5.8)will produce the correct solution in any case that may be 
expanded as a superposition of plane waves. 

We note that the constant terms in the initial data expression (5.4), (5.5)play 
a very important role in the formalism and cannot be ignored. 
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7. Relation to the Hamiltonian Viewpoint 

The evolution operators in (5.8)are close symbolic analogues of the usual evo- 
lution operator in expression (2.3), as we will discuss in some detail. 

A Taylor series expansion in t of a function may be expressed symbolically as 
an exponential as follows 



That is the quantity in the exponent generates a time shift. In standard quan- 
tum mechanics the time derivative is related to the energy operator by ^ = —iE. 
To define the dynamics the energy operator in terms of t is identified with the 
Hamiltonian operator H in terms of x. That is a Schroedinger equation is postu- 
lated, with an operator equivalence 



Then for a wave function the Taylor series expression reads 



which is the standard form for the time evolution operator, as in (2.3). 

We can use the same procedure to obtain our result (5.8)- but only symbolically 
and up to a normalization and constant terms. We begin with an expression for a 
double Taylor series expansion of a function of u and v in analogy with (7.1). 
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(7.1) 



n=0 




(7.3) 



i;{u, v) = i[e("-^°)^^(ii, ^;o) + e(")^V^(Mo, v)] (7.4) 

In analogy with the Schroedinger equation we define the dynamics with the 
Klein Gordon equation, which in terms of u and v is, 

ip{u,v) — — —ip{u,v) (7.5) 



dudv 4 

- which solutions of the Dirac equation must obey. In analogy with (7.2)this 
allows us to identify the v derivative operator with the following operator in terms 
of u, 

d___rrl ,^_m^ 

dv~ A^dJ - 4'" ^^-^^ 

Substituting this into (7.4)we obtain (5.8)- if we do not concern ourselves about 

the normalization factor and the constant terms, which the symbolic derivation 
does not seem to explain. 

In conclusion our main result (5. 8) presents the evolution of a spinor field in 
terms of two analogues of the Hamiltonian operator, and the dynamics of the spinor 
field is remarkably similar to that of the scalar field. 
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